Spin injection across magnetic/non-magnetic interfaces with finite magnetic layers 
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We have reconsidered the problem of spin injection across ferromagnet/non-magnetic- 
semiconductor (FM/NMS) and dilute-magnetic-semiconductor/non-magnetic-semiconductor inter- 
faces, for structures with finite magnetic layers (FM or DMS) . By using appropriate physical bound- 
ary conditions, we find expressions for the resistances of these structures which are in general different 
from previous results in the literature. When the magnetoresistance of the contacts is negligible, we 
find that the spin-accumulation effect alone cannot account for the d dependence observed in recent 
magnetoresistance data. In a limited parameter range, our formulas predict a strong d dependence 
arising from the magnetic contacts in systems where their magnetoresistances are sizable. 
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I. INTRODUCTION 

Spin injection across interfaces is one of the cru- 
cial ingredients for the successful implementation 
of novel spintronic devices || H ||. For in- 
stance, the spin-transistor proposal Q relies on sub- 
jecting injected spin-polarized electrons to a con- 
trollable spin precession between the ferromagnetic 
source and drain. Hybrid ferromagnetic/semiconductor 
and dilute-magnetic-semicondutor/non-magnetic semi- 
conductor junctions constitute basic interfaces in which 
to investigate spin-polarized transport. However, the 
efficiency of spin-injection through ideal FM/NMS in- 
terfaces turns out to be disappointingly small due to 
the large conductivity mismatch |f| between the metallic 
ferromagnet and the semiconductor. The use of spin- 
selective interfaces can significantly enhance injection ef- 
ficiencies . This can be accomplished by inserting spin- 
dependent tunnel barriers between the FM and the NMS 
layers 0. 

Particularly promising is spin injection from a dilute- 
magnetic semiconductor (DMS) into a non-magnetic 
semiconductor (NMS). Novel DMS/NMS junctions (i) 
minimize the conductivity mismatch problem and (ii) 
naturally incorporate spin dependency in the transmis- 
sion process As recently demonstrated, substantial 
spin injection can be achieved in these Mn-based het- 
erostructures • More recently, a novel large magnetore- 
sistance effect has been observed in a DMS /NMS /DMS 
geometries 10]. An interesting observation of Ref. 0] 
is the dependence of the magnetoresistance effect on the 
thickness of the DMS layers: the MR effect doubles with 
the DMS thickness. 

Available formulas describing GMR-type effects in 
magnetic/non-magnetic junctions assume ID geometries 
with semi- infinite magnetic layers |ll|. In this case, 



the expressions for the electrochemical potentials for the 
spin-up and the spin-down electrons contain only decay- 
ing exponentials in the magnetic regions |l2]. The re- 
sistances in the parallel and antiparallel configurations 
are calculated assuming that the device extends for a 
spin-flip length into the magnetic contacts. In such an 
approach it is not clear how to properly take into account 
the voltage drop across the sample. 
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FIG. 1: Color online. Schematic of the magnetic/non- 
magnetic structure investigated here (a) and the correspond- 
ing spin-resolved electrochemical potentials (b) versus the po- 
sition across the layers. Note that we consider magnetic lay- 
ers (DMS) of finite widths d symmetrically placed around the 
non-magnetic layer (NMS). The spin-splitting of the electro- 
chemical potentials across the magnetic/non-magnetic inter- 
faces is larger for a higher "conductivity mismatch" between 
these layers. 

In this work we consider DMS /NMS /DMS and 
FM/NMS/FM ID systems with finite magnetic layers, 
Fig. 1. As detailed below, we describe diffusive trans- 
port in these structures via the usual diffusion theory 
of van Son et al. with proper physical boundary 
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conditions between the several magnetic/non-magnetic 
interfaces in the system. The formulas for the magne- 
toresistance obtained within our treatment are in gen- 
eral different from the ones previously obtained in litera- 
ture, even in the experimentally relevant regime d 3> Xd 
and xq < A at; A_d,Aat denote the spin- flip lengths in the 
DMS and NMS layers, respectively, d is the length of each 
DMS layer and 2xq is the length of the NMS layer. How- 
ever, in the case of a FM/NMS/FM system where the 
conductivity mismatch is large, this difference is not im- 
portant since the correction we find is very small in this 
case. For DMS /NMS /DMS systems, on the other hand, 
the conductivity mismatch is small and the magnetore- 
sistance effect can be large. Our formulas significantly 
deviate from the earlier ones in this case, especially in 
the regime where Xd, which can be magnetic-field de- 
pendent, is comparable to d. 



II. SYSTEM AND APPROACH 

Let us denote by 1,2, and 3 the regions correspond- 
ing to the DMS (or FM), NMS, DMS (or FM) layers, 
respectively. Note that in Fig. 1 x = corresponds to 
the center of the structure (NMS layer), x = ±Xq corre- 
spond to the DMS/NMS interfaces and x = ±(xq + d) 
to the metal contact/DMS interfaces. Consider first the 
DMS /NMS /DMS case or, equivalently, a parallel config- 
uration FM/NMS/FM. Our task is to find the resistance 



J 



of the structure R = V/j, where j = j'j + jj_ is the total 
current through the structure and V is an applied bias 
(+V/2 is the potential of the left metal contact and — V/2 
is the potential of the right contact). 

Diffusive transport equations. Our starting point is 
the diffusive transport approach of Ref. 0, with the ba- 
sic equations relating the spin-dependent electrochemical 
potentials and current densities 



(i) 



which hold in all three regions 1, 2, and 3 (see Fig. 1) of 
our system; cr-^j denotes the conductivity of the spin up 
(or down) electrons in the respective layer. The difference 
in the spin-dependent electrochemical potential is |13| 



Ml - Ml = X' 



d 2 (M - Ml) 
dx 2 



(2) 



where A is the spin- flip length (Xd or A at) in the corre- 
sponding layer. As usual, to solve the problem we first 
assume a general solution in each of the three regions 
and then match them at the many interfaces via proper 
boundary conditions. Below we outline this procedure, 
clearly stating the physical boundary conditions we use 
and the simplifying symmetries of the problem. 

General solution for /if 4. Consistent with Eq. 
we use the general expression for the electrochemical po- 
tentials in region 1 (DMS) 



Ml 
Ml 



(Axx + Bx) 



+ (CiX D exp(x/ X D ) + DiX D exp(-x/ X D )) 



1/of 



(3) 
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Here A = A^ is the spin flip length in the DMS layer. 
Note the presence of the exponentially-growing term in 
Eq. J3J, since we consider a finite DMS layer. The cor- 
responding expressions for the electrochemical potentials 
in the regions 2 and 3 are obtained from Eq. (J3J) by the 
replacements A\ — > A 2 ^, B\ — * -62.3, C\ —* C23, and 
D\ — > D-2,3', in addition we should make er;j?i — > a = ctn /2 
and Ad — > Aat in the NMS region 2. The chosen form 
of the solution guarantees conservation of the total cur- 
rent everywhere in the sample. From the symmetry 
of the problem in the parallel configuration considered 
here, it follows that the electrochemical potentials /if 
and fii should be odd functions of the coordinate x, see 
Fig. 1(b). In terms of the coefficients this means that 
A 3 = A U B 3 = -B l7 C 3 = -Di,D 3 = -Ci,B 2 = 0, and 
C2 = —D-2- Hence we can consider only the boundary 
conditions at the left metal contact/DMS interface and 
left DMS/NMS interface. 

Boundary conditions. At the metal contact/DMS in- 
terface we use the continuity of the total current and the 
equality of the electrochemical potentials: /Xf = Hi = 



—eV/2. The condition /if = /tj is consistent with assum- 
ing a metal contact with infinite conductivity. Thus we 
obtain the following set of equations: 



C*i + D 1 evp(2(x + d)/X D ) = 0, 



-A 1 (x + d) + B 1 = V/2, 



and 



J 



-A^a 



(4) 



(5) 



(6) 



At the DMS/NMS interface we use the continuity of 
the electrochemical potentials, /if(— Xq) = /x-f(— %q), 



Ml(- 



"), the conservation of the total cur- 



rent j = jf + j± , and conservation of the jj component 
across the interface. Then we have the following addi- 
tional set of equations 
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A x xq + B l 



X D /af 



r D ) (Ci exp(-^ /A£i) + Di gxp(xo/Ad)) = -A 2 a:o T 2C 2 Aat sinh(:Eo/AAr)/cr, (7) 



-Aid? - (Ci cxp(-a;o/A_ D ) - L>i cxp(a; /A D )) = -A 2 a - 2C 2 cosh(a;o/Ajv), 



(8) 



and 



here. 



3 = -2aA 2 . 



(9) 



Hence we have seven equations, Eqs. Q-®, for seven 
unknown quantities Ai, A2, Bx, C%, C2, D\, and j. By 
solving these equations we can determine the resistance 
of the structure. This we do next. 



III. RESULTS 

For a FM/NMS/FM and a DMS/NMS/DMS system 
( "parallel" configuration) , we find for R p = V/j 



2xo 

rrN 



2d 

~~D 



a 



2/3 2 



(l-^)^coth(^) + ^coth(^)' 



(10) 



where ft = (af - of)/{af + of), <r N = 2<r is the total 
conductivity in the NMS layer and a D — erf + erf is the 
conductivity in the DMS layer; d is the length of each 
DMS layer and 2x is the length of the NMS layer. The 
term proportional to the magnetic-layer width d in Eq. 
(I10JI is not present in earlier formulas in the literature. 
Note that the term proportional to /3 2 in Eq. (|1L)|) differs 
from that in Eq. (1) of Ref. ^(j because it is d dependent 
and also because of its distinctive xq dependence 14]. 

For the antiparallel configuration (opposite FM lay- 
ers), we can similarly obtain an expression for the resis- 
tance R ap - Interestingly, R ap can be obtained from R p 
in Ea. lfTUI) by changing tanh(xo/AAr) to coth(x /AAr) 



Ran — 



2x 



2d 

~D 



a 



2[3 2 



(l-/^coth(^) + ^tanh(fM 



(11) 



In obtaining Eq. (|llfl we have also assumed that the 
spin- up and the spin-down conductivities of region 3 are 
erf and erf, respectively. We emphasize that the term 
proportional to /3 2 [Eqs. lfTU|) . (fTTfl ] does not reduce to 
that in Ref. for d — > 00 since the boundary conditions 
used here are different. 

Magnetoresistances. For metallic FM/NMS/FM sys- 
tems we define the resistance change as usual, i.e., Ai? = 



However, in the physically relevant limit <C d 
1. 
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FIG. 2: Color online. AR/R P as a function of /3 for 
FM/NMS/FM structures with finite FM layers d. The solid 
and dashed curves correspond to the formula derived in the 
present paper [Eqs. (HOP and llll l. while the diamond curve 
is obtained from Eq. (8) of Ref. [3 Here we use xo = 1/im, 
Ad = 10 nm, Ajv — > 00, d — 100 /im, and a D — lOOu^. 
The three curves do not differ significantly for the parameters 
used. 



and xo <C Xn we find from Eqs. (|10H and (|11|) 



AR 

R„ 



D\2 



(3 2 {a N \ D /xvo u ) 



(l-^Kl-^ + ^T^) 2 ^ 



(12) 

In deriving the above equation we have approximated R p 
in the denominator by its dominant contribution 2xq j a 
from the central NMS layer. Again, our Eq. (|12fl is 
different from the corresponding one [Eq. (8)] in Ref. 5]. 
However, as can be seen from comparison of Eq. (|12J) and 
Eq. (8) of Ref. |5| , this deviation is not really important 
for metallic ferromagnets. Indeed, AR/R p which follows 
from our new equations practically coincides with the 
quantity in Eq. (8) of Ref. ^ for realistic values of the 
parameters when a N \d/v d \n *C 1, <t n \d/o' D xo "C 1, 
and for /3's not anomalously close to unity. 

In DMS/NMS/DMS systems, on the other hand, only 
the parallel configuration is possible. In this case, we de- 
fine the relevant resistance change in the system with 
respect to the zero magnetic-field resistance Rq, i.e., 



Rap — Rp- The full expression is too lengthy to be shown Ai? — R p — Rq. We find 
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AiT 

Rq 



(l-/F)£coth(£) + &coth(&) 



r 



d(l/a D (H)-l/a D (0)) 



+ 



(13) 



Note that R is simply the combined NMS/DMS layer 
resistances 2xo/a N + 2d/a D , since no spin accumulation 
occurs for zero magnetic field. In Eq. I|13ll we have indi- 
cated explicitly the additional magnetoresistance arising 
from the DMS (or FM) regions through a D {H). 

Note that AR D /Ro goes to zero for either d — > oo 
or xq — > oo. The exact expression for AR/R p obtained 
from Eqs. i|10|) and i|ll|) [not the approximate one in 
Eq. I)12[l] also vanishes in these limits. Note that this 
vanishing of the resistance change is symmetrical in d 
and xo as expected. Earlier formulas in the literature 
do not show this symmetry since they assume different 
boundary conditions. 

In contrast to metallic FM/NMS/FM systems, 
DMS/NMS/DMS structures can have conductivities of 
the DMS and NMS regions which are comparable. Be- 
sides, A_d can be of the same order as d in some experi- 
ments. Experimental results in this case, e.g., the d de- 
pendence of the magnetoresistance effect, should be inter- 
preted with the use of Ea. i|13|) . As we mentioned earlier, 
the magnetoresistance effect observed in Ref. [l0| dou- 
bles with the DMS thickness. Let us find out whether our 
ID formula can explain this observed thickness depen- 
dence. Note that the d-dependence in Eq. (|T3)l appears 
through the coth(rf/ X D ) term (spin accumulation effect). 
Since a D = aP + aP can in general be magnetic-field 



dependent, there is also a possible d-dependence arising 
from the second term which describes the resistances of 
the DMS leads. 

In order to investigate in more detail the d-dependence 
in AR D /Rq, we need to know the magnetic- field depen- 
dences of the quantities entering AR d /Rq [Eq. 1|13|) ]. 
particularly the magnetic field dependence of the spin- 
dependent conductivities aP ^. Unfortunately, these de- 
pendences are not known. However, we can still gain 
some insight into the problem by neglecting the mag- 
netic field-dependence of the conductivities. In this case, 
the second term in Ea. l|13|) is zero and we can investi- 
gate the expected d-dependence arising from the spin- 
accumulation effect alone. Note that the quantity f3 is 
null in the absence of a magnetic field and increases as 
a function of it. Hence, we plot AR D / Rq as a function 
of degree of spin-polarization (3 for different values of d. 
Figures 3(a), 3(b) show plots of AR D /R vs /3 for dif- 
ferent parameters. We have chosen the values of d, xq 
and <j d / a N = 1/3 corresponding to the real experimen- 
tal situation of the Ref. 01- From these plots we can 
see that even though the absolute value of the magne- 
toresistance is comparable to the experimental value, the 
d-dependence is not as pronounced as in the experiment 
of Ref. . It is indeed true that the magnetoresistance 



effect is larger for increasing d. However, this behavior 
only appears at values of Ad comparable to the thickness 
d and even for \d = 0.5 /im > d, Fig. 3(b), the effect 
is way too small to explain the experimental data. We 
conclude that the spin accumulation effect by itself does 
not explain - with realistic parameters - the observed 
d-dependence of the data in Ref. [Tc| . 





1 




.li 




.le-l-j 


AR D 


.le-2i 


Ro 


.le-3 n 




le-05i 




le-06 : 




1.1 




.li 




.le-l-j 


AR D 


.le-2i 




.le-3 ; 




le-05i 




le-06 : 



(a) °i>/ 






= 10//m 




= 2 jUm 


f K 


= 0.05 jurn 




0.2 0.4 0.6 0.8 

P 



FIG. 3: Color online. AR D /Ro vs (5 for two distinct DMS 
thicknesses and spin-flip lengths. Only the spin accumulation 
part of the magnetoresistance, i.e., the first term in Eq. 1131 . 
is considered here. For d 3> A_d the magnetoresistance does 
not show any dependence on d (a). On the other hand, for Ad 
comparable to d the magnetoresistance (b) increases slightly 
with d; this effect is too small to explain the d-dependence of 
the magnetoresistance in Ref. [To|| . 

Let us mention here the possibility to explain the ex- 
periment by the magnetoresistance of the leads, see the 
second term in Eq. (|13f) . Control experiments on a layer 
of the DMS injector material patterned into a Hall ge- 
ometry showed ~ 1% negative MR [To| . However, the 
transport problem we deal with here can be shown to re- 
duce to that of a ID system |hJ. Note that in the case 
when the magnetic field dependence of a D is well pro- 
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nounced and for realistic values io > ^ > an d for 
ft ~ 1 the first term in Eq. (|13|l is estimated as \n/a D 
while the second term is of the order of d/a D . Thus in 
the case when the magnetoresistance of the leads is not 
zero the d-dependence associated with it can be more 
pronounced than the one associated with the spin accu- 
mulation effect. 

Some limiting cases. From the general equations we 
derived above, we can obtain expressions for some inter- 
esting limiting cases. If a® — > 0, i.e. ft — ► 1, we find from 
Eq. (HDJ) 



„ 2xn 2d Aw , , , , , , 
Rp = — w H fj + 2^ tanh(x /\ N ). (14) 



If, in addition, Ajy — > oo, then 



2x 2d 



2xo 



Hence the resistance of the central part (NMS) doubles 
as expected. Note that the resistance of the DMS part 
does not necessarily double. 

Current spin polarization. We can also easily deter- 
mine some other relevant quantities in our system. For 
example, the degree of spin polarization a(x) of the cur- 
rent at the DMS/NMS (or FM/NMS) interface (x = x ) 
is equal to 



a(x ) = 



it +h i + 



ft 



("1 — ^ 2% l A « tanh(a:o/A]v) 
"V 1 P I \ D tanhfd/An) 



, (15) 



for parallel alignment. This result coincides exactly with 
Eq. (7) of Ref. [| in the limit xq <C A at, Xd <C d. For 
metallic ferromagnetic contacts we can have antiparal- 
lel alignment as well. Note that a(x) is an even (odd) 
function of x for parallel (antiparallel) alignment. 

For completeness, we give below some other expres- 
sions for a(x). At x — ±(xo + d) (even function) we have 



a(x = — xq— d) — ft- 



ft 



cosh(d/A D )[l + g^ ^^iS^) ] 

(16) 

Note that for d 3> Xd ot(x = —xq — d) approaches the 
maximal possible value ft. In addition 



a(x = 0) = a(x = —xq) 



1 



cosh(a;o/Ajv) 



(17) 



Thus in the middle of the sample the spin current po- 
larization is exponentially suppressed. In the case of an- 
tiparallel configuration the corresponding values of a at 
x = —xq and x — —xq — d can be obtained from formulas 
(|T3|l and (JTJJ, respectively, by replacing tanh(xo/AAr) for 
coth(a;o/AAr). 

In conclusion, we have derived the expressions for 
the resistances of FM/NMS/FM (or DMS/NMS/DMS) 
structures using physical boundary conditions at the cor- 
responding interfaces. We have considered magnetic lay- 
ers (DMS or FM) of finite length d. When the magne- 
toresistance of the contacts is negligible (i.e., when or> is 
not magnetic field dependent), we find that the d depen- 
dence of the spin-accumulation contribution to the mag- 
netoresistance is not enough to explain the experimental 
results of Ref. ^(j. On the other hand, in systems where 
ctd is magnetic field dependent our formulas allow for 
a sizable d dependence arising from the magnetoresis- 
tance of contacts. Interestingly, for realistic parameters 
xq d Xd and with ft ~ 1 the magnetoresistance of 
the contacts is larger than that from the spin accumula- 
tion effect. 

This work was supported by the Swiss NSF, NCCR 
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